Include observed data in Figure 7 or similar figure. Show unnormalised. Also quantify in WT and OE the number and fraction of transcripts that are localized. Demonstrate that particles are equivalent between WT and the OE mutant.
To investigate: Inhomogeneous production due to patchy expression of the gal-4 driver. Inhomogeneous transport through ring canals. Waiting of mRNA to form into complexes at the nuclear pore complex, such as for assembly with dynein or other proteins. (Model comparison between some of these are the ‘normal’ model).
To illustrate these ideas, we plot the RNA distributions for WT and OE based on observed data.
The same thing but without the normalisation:
These plots of the observed data suggest that the fraction of complexes localized is different between wild type and overexpressor. The number of transcripts localized also appeasrs lower in the overexpression mutant, but in my opinion this is due to the age/size of the egg chambers considered. The egg chambers quantified are generally smaller. (Possible bias here as there are too many particles in the larger overexpression datasets so that MATLAB runs out of memory in running the FISH-quant code to localize the transcripts.)
We fit a linear model to the data on total RNA counts in each egg chamber for each phenotype. If \[z = \sum_{i=1}^{16} y_i\] is the total RNA count in each egg chamber, then \[\frac{\text{d} z}{\text{d} t} = 15a.\] Therefore \[ z = z_0 + 15at.\]
Based on the total RNA counts in each egg chamber, we can estimate production rate \(a\) directly for WT and OE data and give a ratio of how much bigger the production for the OE is.
## [1] 9.714501
## [1] 24.57184
## [1] 2.529398
This neglects effects from assembly in the oocyte, which we ought to account for.
## [1] 17.74595
## [1] 32.81549
## [1] 1.849182
Based on this it seems reasonable to keep using \(2a\) as the production in the overexpressor.
## [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13]
## [1,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [2,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [3,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [4,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [5,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [6,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [7,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [8,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [9,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [,14] [,15] [,16]
## [1,] 2 2 2
## [2,] 2 2 2
## [3,] 2 2 2
## [4,] 2 2 2
## [5,] 2 2 2
## [6,] 2 2 2
## [7,] 2 2 2
## [8,] 2 2 2
## [9,] 2 2 2
## [,1] [,2] [,3]
## [1,] 1 1 1
## [2,] 1 1 1
## [3,] 1 1 1
## [4,] 1 1 1
## [5,] 1 1 1
## [6,] 1 1 1
## [7,] 1 1 1
## [8,] 1 1 1
## [9,] 1 1 1
##
## Computed from 4000 by 9 log-likelihood matrix
##
## Estimate SE
## elpd_loo -6719.1 5505.6
## p_loo 766.9 624.1
## looic 13438.3 11011.2
## ------
## Monte Carlo SE of elpd_loo is NA.
##
## Pareto k diagnostic values:
## Count Pct. Min. n_eff
## (-Inf, 0.5] (good) 1 11.1% 788
## (0.5, 0.7] (ok) 1 11.1% 354
## (0.7, 1] (bad) 1 11.1% 171
## (1, Inf) (very bad) 6 66.7% 1
## See help('pareto-k-diagnostic') for details.
## [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13]
## [1,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [2,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [3,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [4,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [5,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [6,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [7,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [8,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [9,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [,14] [,15] [,16]
## [1,] 2 2 2
## [2,] 2 2 2
## [3,] 2 2 2
## [4,] 2 2 2
## [5,] 2 2 2
## [6,] 2 2 2
## [7,] 2 2 2
## [8,] 2 2 2
## [9,] 2 2 2
## [,1] [,2] [,3]
## [1,] 6 9 1
## [2,] 16 1 1
## [3,] 12 16 1
## [4,] 10 11 1
## [5,] 16 1 1
## [6,] 10 1 1
## [7,] 1 1 1
## [8,] 6 11 16
## [9,] 6 11 1
##
## Computed from 4000 by 9 log-likelihood matrix
##
## Estimate SE
## elpd_loo -6670.7 5497.0
## p_loo 758.8 644.3
## looic 13341.4 10994.0
## ------
## Monte Carlo SE of elpd_loo is NA.
##
## Pareto k diagnostic values:
## Count Pct. Min. n_eff
## (-Inf, 0.5] (good) 1 11.1% 774
## (0.5, 0.7] (ok) 2 22.2% 283
## (0.7, 1] (bad) 1 11.1% 15
## (1, Inf) (very bad) 5 55.6% 1
## See help('pareto-k-diagnostic') for details.
## [1] TRUE TRUE
## [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13]
## [1,] 0 1 1 1 1 1 1 1 1 1 1 1 1
## [2,] 0 1 1 1 1 1 1 1 1 1 1 1 1
## [3,] 0 1 1 1 1 1 1 1 1 1 1 1 1
## [4,] 0 1 1 1 1 1 1 1 1 1 1 1 1
## [5,] 0 1 1 1 1 1 1 1 1 1 1 1 1
## [6,] 0 1 1 1 1 1 1 1 1 1 1 1 1
## [7,] 0 1 1 1 1 1 1 1 1 1 1 1 1
## [8,] 0 1 1 1 1 1 1 1 1 1 1 1 1
## [9,] 0 1 1 1 1 1 1 1 1 1 1 1 1
## [,14] [,15] [,16]
## [1,] 1 1 1
## [2,] 1 1 1
## [3,] 1 1 1
## [4,] 1 1 1
## [5,] 1 1 1
## [6,] 1 1 1
## [7,] 1 1 1
## [8,] 1 1 1
## [9,] 1 1 1
## [,1] [,2] [,3]
## [1,] 1 1 1
## [2,] 1 1 1
## [3,] 1 1 1
## [4,] 1 1 1
## [5,] 1 1 1
## [6,] 1 1 1
## [7,] 1 1 1
## [8,] 1 1 1
## [9,] 1 1 1
##
## Computed from 4000 by 9 log-likelihood matrix
##
## Estimate SE
## elpd_loo -7116.7 5817.1
## p_loo 983.3 766.9
## looic 14233.4 11634.3
## ------
## Monte Carlo SE of elpd_loo is NA.
##
## Pareto k diagnostic values:
## Count Pct. Min. n_eff
## (-Inf, 0.5] (good) 1 11.1% 850
## (0.5, 0.7] (ok) 0 0.0% <NA>
## (0.7, 1] (bad) 1 11.1% 46
## (1, Inf) (very bad) 7 77.8% 1
## See help('pareto-k-diagnostic') for details.
## [1] 2 2
## [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13]
## [1,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [2,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [3,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [4,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [5,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [6,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [7,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [8,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [9,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [,14] [,15] [,16]
## [1,] 2 2 2
## [2,] 2 2 2
## [3,] 2 2 2
## [4,] 2 2 2
## [5,] 2 2 2
## [6,] 2 2 2
## [7,] 2 2 2
## [8,] 2 2 2
## [9,] 2 2 2
## [,1] [,2] [,3]
## [1,] 1 1 1
## [2,] 1 1 1
## [3,] 1 1 1
## [4,] 1 1 1
## [5,] 1 1 1
## [6,] 1 1 1
## [7,] 1 1 1
## [8,] 1 1 1
## [9,] 1 1 1
##
## Computed from 4000 by 9 log-likelihood matrix
##
## Estimate SE
## elpd_loo -6596.8 5534.3
## p_loo 780.0 727.8
## looic 13193.5 11068.5
## ------
## Monte Carlo SE of elpd_loo is NA.
##
## Pareto k diagnostic values:
## Count Pct. Min. n_eff
## (-Inf, 0.5] (good) 3 33.3% 294
## (0.5, 0.7] (ok) 0 0.0% <NA>
## (0.7, 1] (bad) 2 22.2% 199
## (1, Inf) (very bad) 4 44.4% 1
## See help('pareto-k-diagnostic') for details.
## elpd_diff se
## 48.4 19.7
## elpd_diff se
## -397.6 314.1
## elpd_diff se
## -446.0 323.5
## elpd_diff elpd_loo se_elpd_loo p_loo se_p_loo looic
## res_prod4a[[2]] 0.0 -6596.8 5534.3 780.0 727.8 13193.5
## res_block[[2]] -73.9 -6670.7 5497.0 758.8 644.3 13341.4
## res_simple[[2]] -122.4 -6719.1 5505.6 766.9 624.1 13438.3
## res_prod[[2]] -519.9 -7116.7 5817.1 983.3 766.9 14233.4
## se_looic
## res_prod4a[[2]] 11068.5
## res_block[[2]] 10994.0
## res_simple[[2]] 11011.2
## res_prod[[2]] 11634.3
We will consider the following models: 0. Simple (M0) 1. Blocked RCs (M1) 2. Inhomogeneous production (2a,a) (M2) 3. Inhomogeneous production (4a,2a) (M3) 4. Inhomogeneous production (4a,a) (M4) 5. Inhomogeneous production (ka,a) (M5) 6. Blocking and production (M6) 7. Density dependent transport (M7) 8. Density dependent transport and production (M8) 9. Density dependent transport and blocking (M9)
## [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13]
## [1,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [2,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [3,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [4,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [5,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [6,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [7,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [8,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [9,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [,14] [,15] [,16]
## [1,] 2 2 2
## [2,] 2 2 2
## [3,] 2 2 2
## [4,] 2 2 2
## [5,] 2 2 2
## [6,] 2 2 2
## [7,] 2 2 2
## [8,] 2 2 2
## [9,] 2 2 2
## [,1] [,2] [,3]
## [1,] 1 1 1
## [2,] 1 1 1
## [3,] 1 1 1
## [4,] 1 1 1
## [5,] 1 1 1
## [6,] 1 1 1
## [7,] 1 1 1
## [8,] 1 1 1
## [9,] 1 1 1
##
## Computed from 4000 by 9 log-likelihood matrix
##
## Estimate SE
## elpd_loo -6719.1 5505.6
## p_loo 766.9 624.1
## looic 13438.3 11011.2
## ------
## Monte Carlo SE of elpd_loo is NA.
##
## Pareto k diagnostic values:
## Count Pct. Min. n_eff
## (-Inf, 0.5] (good) 1 11.1% 788
## (0.5, 0.7] (ok) 1 11.1% 354
## (0.7, 1] (bad) 1 11.1% 171
## (1, Inf) (very bad) 6 66.7% 1
## See help('pareto-k-diagnostic') for details.
## [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13]
## [1,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [2,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [3,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [4,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [5,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [6,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [7,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [8,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [9,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [,14] [,15] [,16]
## [1,] 2 2 2
## [2,] 2 2 2
## [3,] 2 2 2
## [4,] 2 2 2
## [5,] 2 2 2
## [6,] 2 2 2
## [7,] 2 2 2
## [8,] 2 2 2
## [9,] 2 2 2
## [,1] [,2] [,3]
## [1,] 6 9 1
## [2,] 16 1 1
## [3,] 12 16 1
## [4,] 10 11 1
## [5,] 16 1 1
## [6,] 10 1 1
## [7,] 1 1 1
## [8,] 6 11 16
## [9,] 6 11 1
##
## Computed from 4000 by 9 log-likelihood matrix
##
## Estimate SE
## elpd_loo -6670.7 5497.0
## p_loo 758.8 644.3
## looic 13341.4 10994.0
## ------
## Monte Carlo SE of elpd_loo is NA.
##
## Pareto k diagnostic values:
## Count Pct. Min. n_eff
## (-Inf, 0.5] (good) 1 11.1% 774
## (0.5, 0.7] (ok) 2 22.2% 283
## (0.7, 1] (bad) 1 11.1% 15
## (1, Inf) (very bad) 5 55.6% 1
## See help('pareto-k-diagnostic') for details.
## [1] TRUE TRUE
## [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13]
## [1,] 0 1 1 1 1 1 1 1 1 1 1 1 1
## [2,] 0 1 1 1 1 1 1 1 1 1 1 1 1
## [3,] 0 1 1 1 1 1 1 1 1 1 1 1 1
## [4,] 0 1 1 1 1 1 1 1 1 1 1 1 1
## [5,] 0 1 1 1 1 1 1 1 1 1 1 1 1
## [6,] 0 1 1 1 1 1 1 1 1 1 1 1 1
## [7,] 0 1 1 1 1 1 1 1 1 1 1 1 1
## [8,] 0 1 1 1 1 1 1 1 1 1 1 1 1
## [9,] 0 1 1 1 1 1 1 1 1 1 1 1 1
## [,14] [,15] [,16]
## [1,] 1 1 1
## [2,] 1 1 1
## [3,] 1 1 1
## [4,] 1 1 1
## [5,] 1 1 1
## [6,] 1 1 1
## [7,] 1 1 1
## [8,] 1 1 1
## [9,] 1 1 1
## [,1] [,2] [,3]
## [1,] 1 1 1
## [2,] 1 1 1
## [3,] 1 1 1
## [4,] 1 1 1
## [5,] 1 1 1
## [6,] 1 1 1
## [7,] 1 1 1
## [8,] 1 1 1
## [9,] 1 1 1
##
## Computed from 4000 by 9 log-likelihood matrix
##
## Estimate SE
## elpd_loo -7116.7 5817.1
## p_loo 983.3 766.9
## looic 14233.4 11634.3
## ------
## Monte Carlo SE of elpd_loo is NA.
##
## Pareto k diagnostic values:
## Count Pct. Min. n_eff
## (-Inf, 0.5] (good) 1 11.1% 850
## (0.5, 0.7] (ok) 0 0.0% <NA>
## (0.7, 1] (bad) 1 11.1% 46
## (1, Inf) (very bad) 7 77.8% 1
## See help('pareto-k-diagnostic') for details.
## [1] 2 2
## [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13]
## [1,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [2,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [3,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [4,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [5,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [6,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [7,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [8,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [9,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [,14] [,15] [,16]
## [1,] 2 2 2
## [2,] 2 2 2
## [3,] 2 2 2
## [4,] 2 2 2
## [5,] 2 2 2
## [6,] 2 2 2
## [7,] 2 2 2
## [8,] 2 2 2
## [9,] 2 2 2
## [,1] [,2] [,3]
## [1,] 1 1 1
## [2,] 1 1 1
## [3,] 1 1 1
## [4,] 1 1 1
## [5,] 1 1 1
## [6,] 1 1 1
## [7,] 1 1 1
## [8,] 1 1 1
## [9,] 1 1 1
##
## Computed from 4000 by 9 log-likelihood matrix
##
## Estimate SE
## elpd_loo -6596.8 5534.3
## p_loo 780.0 727.8
## looic 13193.5 11068.5
## ------
## Monte Carlo SE of elpd_loo is NA.
##
## Pareto k diagnostic values:
## Count Pct. Min. n_eff
## (-Inf, 0.5] (good) 3 33.3% 294
## (0.5, 0.7] (ok) 0 0.0% <NA>
## (0.7, 1] (bad) 2 22.2% 199
## (1, Inf) (very bad) 4 44.4% 1
## See help('pareto-k-diagnostic') for details.
## [1] 4 1
## [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13]
## [1,] 0 1 1 1 1 1 1 1 1 1 1 1 1
## [2,] 0 1 1 4 1 4 4 4 1 4 4 4 4
## [3,] 0 1 1 4 1 1 4 4 1 1 4 4 1
## [4,] 0 4 1 4 1 4 1 1 1 1 1 4 4
## [5,] 0 1 1 1 1 1 4 1 1 1 4 1 1
## [6,] 0 1 1 4 1 1 4 4 1 4 4 4 4
## [7,] 0 1 1 4 1 4 4 4 1 4 4 1 4
## [8,] 0 1 1 4 1 1 4 4 1 1 4 4 4
## [9,] 0 1 1 4 1 1 4 4 1 1 4 4 4
## [,14] [,15] [,16]
## [1,] 1 1 1
## [2,] 4 4 4
## [3,] 1 4 4
## [4,] 1 4 4
## [5,] 1 1 1
## [6,] 1 4 4
## [7,] 4 4 4
## [8,] 1 4 4
## [9,] 1 4 4
## [,1] [,2] [,3]
## [1,] 1 1 1
## [2,] 1 1 1
## [3,] 1 1 1
## [4,] 1 1 1
## [5,] 1 1 1
## [6,] 1 1 1
## [7,] 1 1 1
## [8,] 1 1 1
## [9,] 1 1 1
##
## Computed from 4000 by 9 log-likelihood matrix
##
## Estimate SE
## elpd_loo -6932.1 5811.5
## p_loo 782.0 680.8
## looic 13864.1 11623.1
## ------
## Monte Carlo SE of elpd_loo is NA.
##
## Pareto k diagnostic values:
## Count Pct. Min. n_eff
## (-Inf, 0.5] (good) 2 22.2% 1411
## (0.5, 0.7] (ok) 1 11.1% 388
## (0.7, 1] (bad) 3 33.3% 17
## (1, Inf) (very bad) 3 33.3% 1
## See help('pareto-k-diagnostic') for details.
## [1] "fitting difference between producers in OE and WT from data"
## [1] -1 1
## [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13]
## [1,] 0 1 1 1 1 1 1 1 1 1 1 1 1
## [2,] 0 1 1 -1 1 -1 -1 -1 1 -1 -1 -1 -1
## [3,] 0 1 1 -1 1 1 -1 -1 1 1 -1 -1 1
## [4,] 0 -1 1 -1 1 -1 1 1 1 1 1 -1 -1
## [5,] 0 1 1 1 1 1 -1 1 1 1 -1 1 1
## [6,] 0 1 1 -1 1 1 -1 -1 1 -1 -1 -1 -1
## [7,] 0 1 1 -1 1 -1 -1 -1 1 -1 -1 1 -1
## [8,] 0 1 1 -1 1 1 -1 -1 1 1 -1 -1 -1
## [9,] 0 1 1 -1 1 1 -1 -1 1 1 -1 -1 -1
## [,14] [,15] [,16]
## [1,] 1 1 1
## [2,] -1 -1 -1
## [3,] 1 -1 -1
## [4,] 1 -1 -1
## [5,] 1 1 1
## [6,] 1 -1 -1
## [7,] -1 -1 -1
## [8,] 1 -1 -1
## [9,] 1 -1 -1
## [,1] [,2] [,3]
## [1,] 1 1 1
## [2,] 1 1 1
## [3,] 1 1 1
## [4,] 1 1 1
## [5,] 1 1 1
## [6,] 1 1 1
## [7,] 1 1 1
## [8,] 1 1 1
## [9,] 1 1 1
##
## Computed from 4000 by 9 log-likelihood matrix
##
## Estimate SE
## elpd_loo -2711.2 1114.4
## p_loo 224.4 193.3
## looic 5422.5 2228.7
## ------
## Monte Carlo SE of elpd_loo is NA.
##
## Pareto k diagnostic values:
## Count Pct. Min. n_eff
## (-Inf, 0.5] (good) 0 0.0% <NA>
## (0.5, 0.7] (ok) 8 88.9% 147
## (0.7, 1] (bad) 0 0.0% <NA>
## (1, Inf) (very bad) 1 11.1% 1
## See help('pareto-k-diagnostic') for details.
## [1] 4 1
## [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13]
## [1,] 0 1 1 1 1 1 1 1 1 1 1 1 1
## [2,] 0 1 1 4 1 4 4 4 1 4 4 4 4
## [3,] 0 1 1 4 1 1 4 4 1 1 4 4 1
## [4,] 0 4 1 4 1 4 1 1 1 1 1 4 4
## [5,] 0 1 1 1 1 1 4 1 1 1 4 1 1
## [6,] 0 1 1 4 1 1 4 4 1 4 4 4 4
## [7,] 0 1 1 4 1 4 4 4 1 4 4 1 4
## [8,] 0 1 1 4 1 1 4 4 1 1 4 4 4
## [9,] 0 1 1 4 1 1 4 4 1 1 4 4 4
## [,14] [,15] [,16]
## [1,] 1 1 1
## [2,] 4 4 4
## [3,] 1 4 4
## [4,] 1 4 4
## [5,] 1 1 1
## [6,] 1 4 4
## [7,] 4 4 4
## [8,] 1 4 4
## [9,] 1 4 4
## [,1] [,2] [,3]
## [1,] 6 9 1
## [2,] 16 1 1
## [3,] 12 16 1
## [4,] 10 11 1
## [5,] 16 1 1
## [6,] 10 1 1
## [7,] 1 1 1
## [8,] 6 11 16
## [9,] 6 11 1
##
## Computed from 4000 by 9 log-likelihood matrix
##
## Estimate SE
## elpd_loo -6938.2 5818.3
## p_loo 748.3 652.1
## looic 13876.4 11636.6
## ------
## Monte Carlo SE of elpd_loo is NA.
##
## Pareto k diagnostic values:
## Count Pct. Min. n_eff
## (-Inf, 0.5] (good) 2 22.2% 1877
## (0.5, 0.7] (ok) 3 33.3% 143
## (0.7, 1] (bad) 1 11.1% 15
## (1, Inf) (very bad) 3 33.3% 1
## See help('pareto-k-diagnostic') for details.
## [1] 2 2
## [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13]
## [1,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [2,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [3,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [4,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [5,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [6,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [7,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [8,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [9,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [,14] [,15] [,16]
## [1,] 2 2 2
## [2,] 2 2 2
## [3,] 2 2 2
## [4,] 2 2 2
## [5,] 2 2 2
## [6,] 2 2 2
## [7,] 2 2 2
## [8,] 2 2 2
## [9,] 2 2 2
## [,1] [,2] [,3]
## [1,] 1 1 1
## [2,] 1 1 1
## [3,] 1 1 1
## [4,] 1 1 1
## [5,] 1 1 1
## [6,] 1 1 1
## [7,] 1 1 1
## [8,] 1 1 1
## [9,] 1 1 1
##
## Computed from 4000 by 9 log-likelihood matrix
##
## Estimate SE
## elpd_loo -6499.5 5141.6
## p_loo 937.4 698.8
## looic 12999.0 10283.2
## ------
## Monte Carlo SE of elpd_loo is NA.
##
## Pareto k diagnostic values:
## Count Pct. Min. n_eff
## (-Inf, 0.5] (good) 0 0.0% <NA>
## (0.5, 0.7] (ok) 1 11.1% 152
## (0.7, 1] (bad) 1 11.1% 38
## (1, Inf) (very bad) 7 77.8% 1
## See help('pareto-k-diagnostic') for details.
## [1] 4 1
## [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13]
## [1,] 0 1 1 1 1 1 1 1 1 1 1 1 1
## [2,] 0 1 1 4 1 4 4 4 1 4 4 4 4
## [3,] 0 1 1 4 1 1 4 4 1 1 4 4 1
## [4,] 0 4 1 4 1 4 1 1 1 1 1 4 4
## [5,] 0 1 1 1 1 1 4 1 1 1 4 1 1
## [6,] 0 1 1 4 1 1 4 4 1 4 4 4 4
## [7,] 0 1 1 4 1 4 4 4 1 4 4 1 4
## [8,] 0 1 1 4 1 1 4 4 1 1 4 4 4
## [9,] 0 1 1 4 1 1 4 4 1 1 4 4 4
## [,14] [,15] [,16]
## [1,] 1 1 1
## [2,] 4 4 4
## [3,] 1 4 4
## [4,] 1 4 4
## [5,] 1 1 1
## [6,] 1 4 4
## [7,] 4 4 4
## [8,] 1 4 4
## [9,] 1 4 4
## [,1] [,2] [,3]
## [1,] 1 1 1
## [2,] 1 1 1
## [3,] 1 1 1
## [4,] 1 1 1
## [5,] 1 1 1
## [6,] 1 1 1
## [7,] 1 1 1
## [8,] 1 1 1
## [9,] 1 1 1
##
## Computed from 4000 by 9 log-likelihood matrix
##
## Estimate SE
## elpd_loo -6606.0 5418.9
## p_loo 862.1 713.6
## looic 13212.0 10837.9
## ------
## Monte Carlo SE of elpd_loo is NA.
##
## Pareto k diagnostic values:
## Count Pct. Min. n_eff
## (-Inf, 0.5] (good) 0 0.0% <NA>
## (0.5, 0.7] (ok) 0 0.0% <NA>
## (0.7, 1] (bad) 4 44.4% 13
## (1, Inf) (very bad) 5 55.6% 1
## See help('pareto-k-diagnostic') for details.
## [1] "using real data \n"
## [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13]
## [1,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [2,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [3,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [4,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [5,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [6,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [7,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [8,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [9,] 0 2 2 2 2 2 2 2 2 2 2 2 2
## [,14] [,15] [,16]
## [1,] 2 2 2
## [2,] 2 2 2
## [3,] 2 2 2
## [4,] 2 2 2
## [5,] 2 2 2
## [6,] 2 2 2
## [7,] 2 2 2
## [8,] 2 2 2
## [9,] 2 2 2
## [,1] [,2] [,3]
## [1,] 6 9 1
## [2,] 16 1 1
## [3,] 12 16 1
## [4,] 10 11 1
## [5,] 16 1 1
## [6,] 10 1 1
## [7,] 1 1 1
## [8,] 6 11 16
## [9,] 6 11 1
## Inference for Stan model: mrna_transport_density_dependent_with_blocking.
## 4 chains, each with iter=2000; warmup=1000; thin=1;
## post-warmup draws per chain=1000, total post-warmup draws=4000.
##
## mean se_mean sd 2.5% 25% 50% 75% 97.5% n_eff Rhat
## a 12.15 0.04 2.05 8.93 10.74 11.88 13.32 17.02 2112 1
## b 14.55 0.16 6.35 4.26 10.04 13.76 18.53 28.96 1673 1
## nu 0.18 0.00 0.05 0.12 0.15 0.17 0.20 0.29 1111 1
## phi 0.38 0.00 0.05 0.29 0.35 0.38 0.41 0.47 2775 1
## sigma 1.30 0.00 0.14 1.05 1.20 1.29 1.39 1.58 2588 1
## beta 0.02 0.00 0.01 0.01 0.01 0.02 0.02 0.03 1635 1
##
## Samples were drawn using NUTS(diag_e) at Sun Oct 14 23:51:21 2018.
## For each parameter, n_eff is a crude measure of effective sample size,
## and Rhat is the potential scale reduction factor on split chains (at
## convergence, Rhat=1).
##
## Computed from 4000 by 9 log-likelihood matrix
##
## Estimate SE
## elpd_loo -6427.9 5157.8
## p_loo 898.7 733.2
## looic 12855.9 10315.7
## ------
## Monte Carlo SE of elpd_loo is NA.
##
## Pareto k diagnostic values:
## Count Pct. Min. n_eff
## (-Inf, 0.5] (good) 0 0.0% <NA>
## (0.5, 0.7] (ok) 1 11.1% 229
## (0.7, 1] (bad) 2 22.2% 27
## (1, Inf) (very bad) 6 66.7% 1
## See help('pareto-k-diagnostic') for details.
## elpd_diff elpd_loo se_elpd_loo p_loo se_p_loo looic
## res_M5[[2]] 0.0 -2711.2 1114.4 224.4 193.3 5422.5
## res_M9[[2]] -3716.7 -6427.9 5157.8 898.7 733.2 12855.9
## res_M7[[2]] -3788.2 -6499.5 5141.6 937.4 698.8 12999.0
## res_M3[[2]] -3885.5 -6596.8 5534.3 780.0 727.8 13193.5
## res_M8[[2]] -3894.8 -6606.0 5418.9 862.1 713.6 13212.0
## res_M1[[2]] -3959.5 -6670.7 5497.0 758.8 644.3 13341.4
## res_M0[[2]] -4007.9 -6719.1 5505.6 766.9 624.1 13438.3
## res_M4[[2]] -4220.8 -6932.1 5811.5 782.0 680.8 13864.1
## res_M6[[2]] -4226.9 -6938.2 5818.3 748.3 652.1 13876.4
## res_M2[[2]] -4405.4 -7116.7 5817.1 983.3 766.9 14233.4
## se_looic
## res_M5[[2]] 2228.7
## res_M9[[2]] 10315.7
## res_M7[[2]] 10283.2
## res_M3[[2]] 11068.5
## res_M8[[2]] 10837.9
## res_M1[[2]] 10994.0
## res_M0[[2]] 11011.2
## res_M4[[2]] 11623.1
## res_M6[[2]] 11636.6
## res_M2[[2]] 11634.3
## elpd_diff elpd_loo se_elpd_loo p_loo se_p_loo looic
## res_M9[[2]] 0.0 -6427.9 5157.8 898.7 733.2 12855.9
## res_M7[[2]] -71.6 -6499.5 5141.6 937.4 698.8 12999.0
## res_M3[[2]] -168.8 -6596.8 5534.3 780.0 727.8 13193.5
## res_M8[[2]] -178.1 -6606.0 5418.9 862.1 713.6 13212.0
## res_M1[[2]] -242.8 -6670.7 5497.0 758.8 644.3 13341.4
## res_M0[[2]] -291.2 -6719.1 5505.6 766.9 624.1 13438.3
## res_M4[[2]] -504.1 -6932.1 5811.5 782.0 680.8 13864.1
## res_M6[[2]] -510.3 -6938.2 5818.3 748.3 652.1 13876.4
## res_M2[[2]] -688.8 -7116.7 5817.1 983.3 766.9 14233.4
## se_looic
## res_M9[[2]] 10315.7
## res_M7[[2]] 10283.2
## res_M3[[2]] 11068.5
## res_M8[[2]] 10837.9
## res_M1[[2]] 10994.0
## res_M0[[2]] 11011.2
## res_M4[[2]] 11623.1
## res_M6[[2]] 11636.6
## res_M2[[2]] 11634.3
TO DO: Run M9 (Done) Rerun M5 to figure out what went on here (or leave this). Run some/all ommitting young egg chambers or investigating timescale model to try to avoid issues due to this